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Abstract 

The charge density and pair correlation function of three interacting electrons confined 
within a two-dimensional disc-like hard wall quantum dot are calculated by full numeri- 
cal diagonalization of the Hamiltonian. The formation of a Wigner-molecule in the form 
of equilateral triangular configuration for electrons is observed as the size of the dot is 
increased. 
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1 Introduction 



Advances in nanostructure technology have allowed the lateral confinement of a two dimensional 
electron gas (2DEG) by means of suitably shaped gate electrodes or by etching techniques [1]. 
These confined systems having a discrete energy spectra are commonly called zero dimensional 
systems or quantum dots [2-9]. The motion of electrons in a quantum dot is quantized in all 
three spatial directions. However, if the quantization in the vertical direction is much stronger 
than the quantization in the in-plane directions, a quantum dot could be treated as disc-like 
2D system where electrons have significant freedom along x- and y- directions. 

Experimentally, the number of electrons confined in a quantum dot could be varied over a 
considerable range by changing the gate voltage applied. Quantum dots containing as few as 2, 
3 or 4 electrons have already been realized and investigated by optical absorption experiments 
[6,9]. In these quantum dots or artificial atoms the Coulomb interaction between the electrons 
is very important for understanding their quantum mechanical properties. Especially correla- 
tions among electrons are crucial since their effects influence the spectral [10-14] and transport 
properties [15-17] of quantum dots. 

So far, various analytical techniques have been devised for handling the electron corre- 
lations in quantum dot systems. One of these techniques include the solution of the many 
particle Schrodinger equation [18]. In this three dimensional approach, interacting electrons in 
non-parabolic quantum dot systems have been investigated using the formalism of Hylleraas 
[19] where interelectron coordinates are built into the wavefunction explicitly. By comparing 
Hylleraas and Hartree-Fock type results, in [18], it is shown how the charge redistributions, 
attributable to correlation interactions, affect ground and excited state energies and electron 
confinement in isolated and coupled two quantum dot systems. One- and two- electron ground 
state energies of a silicon sphere embedded in an amorphous silicon dioxide matrix [20] have 
been calculated as a function of sphere size. The electron-electron interaction and polarization 
effects in that study have been treated by perturbation. 

A great deal of interest has also gone into analytical investigation of correlation effects 
in 2D quantum dot systems. Comparison of energies, pair correlation functions, and particle 
densities of the singlet and triplet ground state of quantum dot helium in a magnetic field has 
been investigated in [14] by using Hartree, Hartree-Fock, and exact diagonalization methods. 
The HF results for the triplet ground state were found to be in very good agreement with 
the exact diagonalization results, proving the importance of the exchange interaction. On the 
other hand, the exact results for the singlet state have disagreements with the HF results. In 
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[14] it has been shown that the disagreement between two results arises due to the electron 
correlations which are neglected in the Hartree-Fock approximation. However, most of the 
work performed on the electron-electron correlations in quantum dots placed in a magnetic 
field [6,8,21], as well as transport experiments [22] and far-infrared spectroscopy [23] has been 
based on a two dimensional quantum dot with parabolic type confinement potential. 

Typically the lateral confinement potential in quantum dots is created by spatially extended 
charge distributions. Therefore, it shows a parabolic characteristic (~ r 2 ) and seems to repre- 
sent fairly well the electrostatically confined electrons. The advantage of using a harmonic type 
of confining potential is the analytical simplicity of the problem since the center of mass and 
relative coordinates decouple. However, in the far-infrared spectroscopy, the radiation field can 
not detect the electron-electron correlations when the confinement is parabolic. This is caused 
by the decoupling of far- infrared radiation with the relative motion of electrons [24] . 

Thus, in order to probe the interaction effects, it has been suggested that the shape of the dot 
should be modified to achieve the coupling of center of mass and relative motions [24]. In [24], 
the heat capacity results have also been presented to suggest that the interaction effects could 
be seen by measuring the thermodynamic properties of the electrons. The magnetization of 
parabolic quantum dots has been also computed and magnetization found to be another probe 
of the interaction effects [25] like the heat capacity. The magnetization of the dots is predicted 
to oscillate with magnetic field because the ground state prefers to be at certain magic values of 
the total angular momentum. This behaviour is explained in a subsequent work [26] as a direct 
consequence of the Pauli principle which enables the electrons to reduce their energy optimally 
only at the magic angular momenta. Recently, vertically coupled quantum dots or artificial 
molecules [27,28] have attracted considerable attention. In [28], a double dot system with three 
spin-polarized electrons was investigated and a sequence of angular-momentum magic numbers 
was found depending on the strength of the interdot tunneling. Besides these, some authors 
have achieved the coupling of center of mass and relative motions of electrons with deviations 
from the exact harmonic confinement [29-31]. With etching techniques or self organized growth 
it is possible to create hard- wall type confinement potential [32]. Since the center of mass 
motion of electrons is coupled with their relative part in this type of confinement potential, the 
effects of correlation influencing energy spectrum, transport, and spectral measurements can 
be studied with the excitation spectra [33]. 

In this work we analyze electron correlation effects in a 2D circular quantum dot with 
hard wall confinement potential. We employ exact diagonalization technique in computing the 
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eigenvalues and corresponding eigenvectors for the ground state numerically. We specialize to 
the case of three spin-polarized electrons with total magnetic quantum number S z = 3/2 and 
total orbital angular momentum M — 0. We investigate charge density and pair correlation 
function as a function of the dot size under these circumstances. With increasing dot size the 
observed formation of Wigner molecule in equilateral triangular configuration of three electrons 
is discussed. In section 2 we give the formalism and details of the calculation. In section 3 and 
4 we report and discuss the results of the numerical calculations. 

2 Model and the method of calculation 

The total orbital angular momentum M = Y^iLi m i, the total spin S and the total magnetic 
quantum number S z = X^Li s i are good quantum numbers for N electrons in the dot due to the 
circular geometry of the dot, and spin independence of the Coulombic interaction. In second 
quantization language, the N electron quantum dot is described by the Hamiltonian 

7~L = ^2 ^K a ^K a K + n ^2 Vk,L,M,N <^K a ^M a LQ<N , (1) 
K ^ K,L,M,N 

where £k = \^n K \m K \ depends only on the n^- th root of the Bessel function of the first kind 
Jm K - In writing eqn. (1), length is measured in units of the dot radius a, and energy in units of 
h 2 1 '(m* a 2 B \ 2 ) , where m* is the electron effective mass determined by the host semiconductor. 
In (1), the Coulomb potential is multiplied by A = a/as which serves as the dimensionless 
coupling constant characterizing the strength of the interaction. Although it is possible to 
perform a perturbative expansion in powers of A to calculate the energy spectrum [34], we 
prefer to diagonalize the Hamiltonian matrix exactly [33,35]. The Coulomb matrix element is 
defined via the one-electron orbitals as 

V k ,l,m,n = J |rf 2 ^xV^(5)^L(^)V(x-x / )^M(x , )^(^), (2) 

where <pa{%) (A=K, L, M, N) are the eigenfunctions of the single free particle Hamiltonian, and 
A is a collective index designating the radial (tia), angular momentum {uia) and spin quantum 
numbers (a a) of the electron. (Pa(x) can be written as 

</? A (x) = <f) nA)rnA {Z)xa A , (3) 

where Xa is the spin wavefunction, and one electron orbital 0„ im (x) has the form 

n , m (x) = -^=- l - ^e ime J H (A; niH |x|) , (4) 
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and the corresponding eigenvalue £k is independent of both spin and and sign of m. As the 
one electron orbitals depend already on the Bessel functions, it is convenient to expand the 
Coulomb potential, V(x — x') = ' * n terms °f the Bessel functions, 

I 00 /*oo 

——= £ / dfce im < , '-«')J w (fcp)J w (V)e-* ( *>-* <) , (5) 
l x ~ x 'l ^=-00^ 

where z> — z< shows the extension of the dot in the longitudinal direction. If the longitudinal 
extension of the dot is a non-negligible fraction of its in-plane size, one is to analyze the 
excitations in this direction, too. However, in the limit of z > — z < — > 0, it is sufficient to probe 
only the transversal plane assuming that the system is in the lowest state for longitudinal 
dynamics. In the following, we follow the latter one and ascribe a small value to z > — z < in the 
calculations. In this approximation the matrix element of the Coulomb potential becomes 

1 1 1 



yn K ,n L ,n M ,n N _ ^ 

\J\m K \+l(kn K ,\rn K \)\ \J\m L \+l{kn L ,\m L \)\ \ J\m M |+1 (^n M ,|m M | , 



^00 rl 

' dk / dp pJ\ mK \{k nK ^ mK \p)J\ mL \{k n ^\ mL \p) 

J 



I ^|mjv|+l (^njv,|mjv 

J\m K -m L \(kp) / ^pV / ^|m M |(^nM,|mAf|PO^|miv|(^n]V,|m. J v|P / )^|m]V-mM|(^P / ) ■ (6) 
J 

One notes that, angular integrations in Vk,l,m,n require m K — m L = m N — m M otherwise 
V k ,l,m,n vanishes. 

For calculating the physical quantities such as energy spectrum, charge density and pair 
correlation function, one needs to compute the matrix elements of the Hamiltonian operator 
between N-electron ground state which can be expressed as 

I*? (7) 
c* 

where |C >= aj^ mi ai ....a\ lN mN aN \0 > is a non-interacting Slater determinant and sum runs 
over all possible configurations of the quantum numbers (n, m, a) satisfying given values of M 
and S z . The expansion coefficients b° c are identified with the eigenvectors of the ground state. 
In this work we are concerned with the quartet ground state \S Z = 3/2, M = > 

\S Z = 3/2, M — >= 4 limiiT ^ m2)T 4 3im3)T |0 >, (8) 

in which all electrons are spin polarized. We will compute all relevant quantities in the quartet 
ground state. 

The charge density of the electrons are defined by 

P (x) = £< ^i^wum^ >, (9) 

a 
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which measures the electron density at a given point x in space. A close inspection of this 
formula reveals that, at a fixed radius, p(x) ~ cos(2M#) so that angular dependence of the 
charge density is determined by the total orbital angular momentum. Hence, for M = angular 
dependence disappears and p(x) assumes only a radial variation dictated by the associated 
Bessel functions. 

Another relevant quantity, the pair correlation function, is a two-point function defined by 

Pc (f,f ) = J2< < N) \<f>U^l^)M^Mm { o N) >, (io) 

cr,cr' 

which is the probability of finding an electron at x given that another one is situated at x'. 

In the next section we perform a numerical computation of the charge density and pair 
correlation function to identify their dependence on the space coordinates as well as the dot 
size A. 

3 Numerical Analysis 

Electron-electron correlations have been computed for a ID hard- wall type dot in [36], and a 
2D parabolic dot in [37]. In the latter the Schroedinger equation is solved for three electrons 
numerically for a fixed dot size. For three electrons with S z = 3/2 and M = 3k (k=0,l,2, ...) 
it is expected that, in the ground state, electrons form an equilateral triangle at some radius 
r determined by the confining parabolic potential. Here we analyze charge density and pair 
correlation function for three electrons in a circular quantum dot for S z = 3/2, M — 0. In 
contrast to the parabolic confinement potential, in the case of hard-wall type confinement one 
cannot find an analitic expression for r . Despite of this, however, for large enough A value, the 
equilateral triangular configuration is expected to occur at some distance in the radial direction. 
This then will be an indication of the Wigner molecule structure which is a general feature of 
quantum dots for large A [36]. 

Depicted in Fig.l is the r and A dependence of the charge density p(x) (normalized to N=3 
for each A) for the state \S Z = 3/2, M = >. The three curves in this figure describe variations 
of the electron distribution as a function of A and r. For small A (e. g. A = 1) electrons have 
non-negligible distribution at the center of the dot though maximum is reached away from the 
origin. This behaviour of p(x) shows that for small A J m=0 (r) is the dominant one especially at 
small r. As A increases, however, the charge density vanishes gradually around the center of the 
dot (e. g. A = 10, 100). It is here that one observes the dominance of J m >o since they vanish 
at the origin by definition. Besides the behaviour of the charge density close to the center of 
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the dot, one observes that higher the A sharper and farther from the center the peaks become. 
Hence, as A increases electrons get shifted towards the periphery (never reaching there due to 
the infinite potential barrier) with a sharper peak showing the most probable radial position 
of the electrons. For example, for A = 100 electrons are most probably distributed along the 
periphery at the radial distance r ~ 0.7. That electrons move towards the periphery with 
increasing A follows from the fact that they try to minimize their electrostatic energy which 
is known to be the dominant component for large A. Behaviour of the charge density gives 
information only on the radial distribution of electrons; therefore, one sould also investigate 
the behaviour of the pair correlation function to obtain the distribution of the electrons in the 
plane of the dot. 

In Fig. 2 the variation of the pair correlation function (normalized to N(N-1)=6) with 9 
and A for r = r' = 0.67 and 9' = is presented. The pair correlation function p c (x, x') gives 
the probability distribution of N — 1 electrons given that one of the N electrons is located at 
x'. For the case of three electrons, the pair correlation function in Fig. 2 shows the angular 
distribution of two electrons along the periphery of the quantum dot. Therefore, the fact that 
pair correlation function has always two distinct peaks and vanishes at 9 = is a restatement 
of the properties of the pair correlation function consistent with the Pauli exclusion principle. 

An important property of the pair correlation function follows from its variation with the 
dot size in units of Bohr radius. One notices that, for small A values peaks are not sharp 
and electrons do not have a well-defined configuration, that is, the pair correlation function 
is not diminished significantly between the two peaks. This small A limit shows nothing but 
the atomic regime of confinement, where dot size is of the order of Bohr radius or smaller, and 
their average kinetic energy exceeds the Coulombic repulsion. It is with the dominant kinetic 
energy of the electrons that they are distributed in the dot without a well-defined configuration. 
It is known that in this limit perturbation theory is a reliable tool to investigate the physical 
parameters of the dot [38], in other words, Hartree-Fock method can be applied just as in the 
few-electron atoms. 

Furthermore, as Fig. 2 shows clearly, when A is increased peaks get sharper. Pair correla- 
tion function gradually assumes two distinct peaks between which there is a strong depletion. 
Therefore, higher the dot size smaller the overlap between the two peaks. Indeed, as A in- 
creases the peaks approach to fixed positions such that their angular seperation is ~ 120°. 
This angular seperation constitutes an equilateral triangular structure as already emphasized 
in other investigations too [37,38]. This geometrical arrangement of the electrons correspond 



7 



to the minimal energy configuration where their kinetic energy is smaller than the electrostatic 
Coulomb energy. The latter is minimized by an equilateral triangular configuration of the elec- 
trons. Approach of the configuration to an equilateral triangular structure for large dot sizes is 
nothing but the well-known Wigner molecule structure which has been shown to exist in other 
types of dots too [36,38]. Therefore, for large dot sizes one observes that the electrons form 
a Wigner molecule, that is, they assume fixed positions in the dot minimizing the dominant 
electrostatic energy. 

For a better understanding of the behaviour of the electron distribution in the dot it may 
be convenient to analyze the pair correlation function in the r — 6 plane in which both radial 
and angular variations become visible. Fig. 3 shows the contour lines of p c (x, x') in r — 6 
plane for 6' — 0, r' — 0.67, and A = 1. In accordance with Fig. 1 and Fig. 2 there is a 
non-negligible correlation between the electrons for small r, also the maxima are reached away 
from the origin, r ~ 0.4. Moreover near the periphery of the dot, the correlation distribution 
is highly suppressed. One also notices that the pair correlation function in between the two 
peaks is not suppressed at all. 

To illustrate the effects of larger dot sizes on Fig. 3 we show in Fig. 4 the pair correlation 
function in r — 9 plane for A = 100. As is seen there are important differences between Fig. 3 
and Fig. 4. First of all, the peaks are pushed towards the periphery of the dot in comparison 
with Fig. 1. Next the peaks are now sharpened and the pair correlation between them are 
reduced significantly. In this sense the difference between Figs. 3 and 4 shows the crytallization 
process of the electrons with increasing Coulombic repulsion among them with growing A. 

4 Conclusion 

In this work we have performed a detailed numerical study of the electron correlation effects 
in a 2D circular quantum dot with hard-wall confinement potential. Our investigations show 
that as the dot size increase gradually configuration of the electrons approach an equilateral 
triangular structure. Hard-wall quantum dots, which may be as interesting as parabolic ones for 
experimental studies, allow for the formation of the Wigner molecule structure of the electrons 
for large dot sizes compared to the Bohr radius. It is known that the optical properties like 
the inelastic light scattering or far infrared absorption [39] and magnetic properties are all 
dependent on the spin of the ground state. Therefore, experimental studies on the quantum 
dots can reveal valuable information on the electron configuration in the dot by concentrating 
on the spin dependent quantities. 
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Figure Captions 



Fig. 1. Variation of the charge density p(x) with radial distance r and the dot size A. 
Fig. 2. Variation of the pair correlation function with 9 and A for r = r' = 0.67 and & = 0. 
Fig. 3. Contour lines showing the pair correlation function on (r — 9) plane for r' = 0.67, 
9' = 0, and A = 1. 

Fig. 4. The same as in Fig. 3 but for A = 100. 



11 



2.5 ,— r 




0.2 0.4 0.6 0.8 1 

r 



Fig.l 



12 




Fig.3 
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Fig.4 
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